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<N Abstract 

We model the observed J ++ mesonic mass spectrum in terms of energies for tightly knotted and 
linked chromoelectric QCD flux tubes. The data is fit with one and two parameter models. We 
predict a possible new state at approximately 1190 MeV and a plethora of new states above 1690 
Q MeV. 
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I. INTRODUCTION 



In 1867 Lord Kelvin suggested [T] that elementary particles can be knots. While his idea of 
knotted fluid vortices in the aether as fundamental objects of nature was revolutionary for his 
time, our present experimental knowledge does not agree with this conjecture. Specifically, 
we now know that knotted fluid vortices are unstable and worse still, the aether does not 
exist. Nevertheless, the idea is attractive for its simplicity of relating fundamental physical 
and mathematical objects, and it should not be discarded out of hand before being tested on 
various other physical systems. One such system where stable knotted configurations may 
exist is quantum chromodynamics (QCD), the subject of the present study. 

Observations from many experiments can be interpreted as signatures of unusual mesonic 
states, i.e., bosonic hadrons that are not pure qq [2]. Such states can be broadly divided into 
the following types: (1) hybrids — bound states of quarks and gluons, like qqG with quan- 
tum numbers J PC = h , 1 h , 1 , 2 h , . . .; (2) exotics — for example, four and six quark 
states, such as qqqq and qqqqqq with quantum numbers J PC = , + ~, 2 + ~, . . .; (3) 
glueballs — states with no valence quarks at all, composed of pointlike or collective glue, e.g., 
string loops a la Nielsen-Olesen [3J, or closed flux tubes. (Even though glueballs do not 
contain valence quarks, there are certainly sea (virtual) quarks within a glueball.) Glueballs 
are among the most studied and least understood classes of particles in QCD [I]. Lattice 
calculations, QCD sum rules, electric flux tube models, and constituent glue models lead to 
a consensus that the lightest valence quark-free state is a glueball with quantum numbers 
J ++ = ++ [5J. On the lattice, usually only a single glueball state below ~2 GeV is con- 
sidered, since all the excitations are expected to be above this energy [Bj, however, a full 
study of topological operators responsible for knots and links on the lattice is computation- 
ally challenging at present and has not yet been carried out. Nevertheless, studying such 
configurations would be interesting and potentially important for a better understanding of 
QCD 0. 

Although glueballs, and the / states they are associated with, are one of the most widely 
discussed problems in hadronic physics, and, while many glueball models have been pro- 
posed [8J, there is still no consensus of what constitutes a glueball beyond its J ++ quantum 
numbers. Here we take an egalitarian approach. Specifically, we model all J ++ mesonic 
states, i.e., all fj and fj states listed by the Particle Data Group (PDG) [2] (/ states, for 
brevity), as knotted or linked chromoelectric QCD flux tubes [10]. Hence we will use the 
term "glueball" loosely as a shorthand for any J ++ state in QCD. 

II. MODEL 

In [11] two of us argued how to generalize various classical ideas from plasma physics to a 
semiclassical model of knotted and linked configurations in QCD. To this end, we first recall 
that Maxwell's equations for an ideal plasma imply that flux lines are locked into the plasma 
flow. This means that if the flux lines are knotted or linked, then the flux line topology 
is conserved as the flow evolves. A classical consequence of conservation of topology is the 
concept of helicity and its conservation in a plasma [12]. Helicity in this context corresponds 
to the degree of Gaussian linking of QCD flux tubes. Keeping only the helicity fixed, while 
allowing other topological changes in flux lines, and minimizing the energy leads to the so- 
called Taylor states [13] in plasma physics. Keeping flux tube topology fixed leads to tight 
knots [H] in QCD. 
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The simplest configuration of linked flux tubes has the form of a Hopf link (denoted 2\) 
in which two unknots are linked together in the simplest way such that the Gaussian linking 
number is one. This is the tight Hopf link, in which fixed diameter d flux tubes, carrying one 
flux quantum each, have the shortest total length. The ratio of the total length of the tubes 
to their diameter d is invariant for all such tight Hopf links. We define the "knot energy" of 
the tight configuration K as its dimensionless length Ik, 

e (K) = l K /d, 

where in this example K = 2f, so then Sq(2\) = 2(2nd)/d = Air. By "knot energy," mathe- 
maticians generally mean some kind of repulsive electrostatic energy such as O'Hara's energy 
[T5] or the Freedman-He-Wang [16] "Mobius energy." Our Eq(K) is called the ropelength of 
the configuration and the configurations which minimize the ropelength are called "tight" or 
"ideal" knots [17]; see also [18]. For the Hopf link and a family of other simple links, where 
the tight configuration is known, see [19]. As expected, for the Hopf link this configuration 
consists of two linked circles passing through each other's centers. 

The Hopf link is the simplest link, but there is an infinite family of topologically different 
link types. Many of these configurations have been tabulated by mathematicians [20j EI] 
and there are different topological invariants (such as linking numbers) which distinguish 
them. However, regardless of which invariants are used to identify the configurations, the 
"frozen-in field" hypothesis implies that tube topology is conserved. Hence, we conjecture 
that the ground states of all systems of flux tubes, with any type of nontrivial linking, are 
the states with the shortest length tubes. 

Similar to linked flux tubes, we also consider self-linked (knotted) flux tubes, where the 
associated tight knot states and dimensionless knot lengths are analogously defined. The 
simplest example of a nontrivially knotted flux tube has the form of the trefoil knot 3i. This 
configuration has knot energy which has been numerically calculated to be £o(3i) ~ 16.3715 
[22] . Note that there are no known analytic forms for the lengths of any tight knots or links 
with nonplanar elements. 

To begin the description of the model, we consider a high energy hadron-hadron collision 
in the process of rehadronization, where there are baryons, mesons and quantized fluxes 
confined to tubes. If the tubes are open, with quarks and anti-quarks at their ends, then 
they are excited baryon or meson states. Our interest is in closed tubes which can be self- 
linked (knotted) or linked with each other. As a key part of our model, we identify all the 
/ states as knotted or linked QCD chromoelectric flux tubes. The topological quantum 
numbers (or knot/link type of the configuration) are what stabilizes the knotted and linked 
configurations, so we assume that non-topological (i.e., unknotted/unlinked) J ++ closed flux 
tube configurations are too unstable to have measurable widths. A configuration with tightly 
knotted or linked flux tubes in the form of the knot or link K will be called f(K). Note 
that topological invariants in QCD typically require instanton and Chern-Simons terms for 
their full description [23], but we will not need these subtleties here. 

As argued above, nontrivial knotting and linking leads to quasi-stable generalized mini- 
mum energy states. This implies the following theorem. 

Theorem 1. For a configuration with a topological charge measured by knotting or linking 
of flux tubes of a constant radius (due to quantized flux, and therefore of constant energy or 
mass per unit length), the generalized minimum energy state is the one that minimizes e{K) 
and therefore ties the knot or link with the shortest tube. 
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Clearly, the minimum energy corresponds to the minimum dimensionless length of tube 
needed to support the topology. Since the length and energy coincide up to a rescaling, the 
proof proceeds trivially by inspection. (Note that this proves a quantum analog of Moffatt's 
1985 conjecture that higher order linking leads to positive lower bounds on configuration 
energy.) We will see below that the approximation of a fixed energy per unit length can be 
improved by an analysis of the effect of field rearrangement within a bent tube. 

We conclude that the quantum case of tight flux tube configurations is much simpler 
than the corresponding classical case where one minimizes energy with a flux constraint. 
However, the quantum case suggests that it may be possible to sum or integrate any large 
number of flux quanta to get the classical result for the generalized minimum energy of a 
Gaussian-linked or higher order topological configuration. (Such a result would complete 
the proof of Moffatt's conjecture for the classical case.) This is a side issue from our main 
purpose here that will be explored elsewhere. 

Let us now proceed with the further description of our phenomenological model, which is 
similar to the model in [11] , but with some minor modifications. Pulling a quark- ant iquark 
(qq) pair apart in the QCD vacuum, a chromoelectric flux tube forms along a path connecting 
the quark and antiquark. If the q is annihilated against the q, or if it is annihilated by a q 
in another qq pair, where the new q is in turn annihilated, etc., to close the path, then a 
closed flux tube containing one flux quantum can form as an unknot or a knot. A flux tube 
following such a curved path could also, for instance, be due to multiple scattering of the 
inital qq pair before its mutual annihilation. If the closed tube is a knot, or if it ends up being 
linked with another closed tube, then such an object has at least one nonzero topological 
quantum number, and this quantity is what tends to stabilize the tight configuration, which 
we either identify with one of the observed / states or use to predict a new state. To allow 
hadronization to run its course, we assume that the typical time scale needed to reach a 
tight configuration is shorter than the typical lifetime of the hadron. 

Since the publication of [11] in 2003, there has been slow continuous physics progress in 
the refinement of the / states data as summarized by the PDG [2], with smaller error bars 
from better statistics and a few new states now listed in their summary tables. The change 
in the / state data that affects our model the most is the PDG's realignment of the mass 
of the a, formerly called the / (600) and now reassigned as the / (500). While there has 
been no new data since 2007, there has been several new reanalyses of compilations of the 
existing data. The 2010 PDG value of the /o(600) mass was reported to be in the range 
400 — 1200 MeV and we previously used 800 ± 400 MeV as our approximation. Now the 
PDG is reporting a mass range for the /o(500) of 400 — 550 MeV. This is a drastic change 
since the range has contracted by a factor of 5 and the central value has dropped by over 
300 MeV. 

Contrary to the rather gradual physics progress, the mathematical knowledge of tight 
knots has changed dramatically from what was used in [11] to what it is today [27]. In 
2003 the lengths of only a handful of tight knots and non-planar links were known, and 
some of those only to an accuracy between 5% and 10%. Now we know the complete 
length spectrum of the first several hundred tight prime knots and non-planar links with an 
accuracy that is assumed to be in the 0.1% — 1.0% range for most of this spectrum. The 
lengths of the composite knots have recently appeared [2H]. We present new computations 
covering the final piece of the puzzle — ropelengths of composite links — for the first time 
below. The current limitation on knot energies needed for the model is due to the fact 
that we are dealing with physical knots and links that can be constricted or distorted (see 
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below), thus increasing the errors on the effective lengths. Even with this caveat, we still can 
advantageously refit the J ++ data by comparing it with high accuracy tight knot and link 
data after adding estimated errors, all of which is collected in Tables [I] through IV below. 
Let us summarize the model assumptions: 

1. There is a one-to-one correspondence between / states and tightly knotted and linked 
chromoelectric flux tubes. 



2. The flux is quantized with one flux quantum per tube. 

3. Knotted and linked flux tubes are stabilized by topological quantum numbers. 

4. The tube diameter is in the ~ 0.1 fm range. (This corresponds to a string tension of 
approximately 400 MeV, which agrees with lattice estimates.) 

5. The quantity J in an fj or fj state is the intrinsic angular momentum of the associated 
knotted solitonic solution of the QCD field equations. 

6. The relaxation to a tight state configuration (via processes where no topology change 
is involved) is faster than its decay rate (via processes with topology change) for an / 

State, i.e., T re i ax <C Tdecay 

One modification from [llj is that we now assume J is the intrinsic angular momentum 
rather than the rotational angular momentum. We do this because the tube diameter is 
now assumed to be smaller and hence the rotational energy level spacing to be larger, ~ 500 
MeV, as opposed to a few MeV for the thicker tubes assumed in [TT]. The other significant 
modification is that we correct the energy due to tube curvature as discussed in the next 
section and include estimated errors due to other physical corrections. 

We begin by identifying ideal knot and link lengths with glueballs and/or predicted 
glueballs, where we include all / states. The lightest candidate is the /o(500), which we 
identify with the shortest knot or link, i.e., the Hopf link 2f; the /o(980) is identified with 
the next shortest knot or link, in this case the trefoil knot 3i, etc. 

Our initial one-parameter fit of the data is shown in Figure [TJ The slope is A tubc = 57 
MeV and \ 2 = 84. The fit is poor mainly because of the constraint imposed by the very 
small error bars on the masses of the ^(1270) and the /i(1285). We will now see how to 
improve the fit when ideal tubes are replaced with physical tubes. Note we are already 
assuming errors of 3% on the knot lengths which anticipates this replacement. 



III. CURVATURE CORRECTIONS 

In the discussion so far we have assumed uniform flux across the cross section of the 
tubes, but flux is not necessarily uniform over the cross section of curved tubes. This leads 
us to define a new energy functional for tubes which we call the "flux tube energy." 

To motivate our definition, we consider the effect of bending on the total energy of a field 
confined to a tube. First recall that the magnetic field of an ideal toroidal solenoid with 
fixed flux falls like 1/p from the symmetry axis. To see this, choose cylindrical coordinates 
(z, p, a) as shown in Fig. [2] and note that symmetry requires the field be in the a direction. 
(For an elementary argument see [29J.) Here we will proceed via a variational argument 
which is a simpler alternative. 
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FIG. 1. Fit using uncorrected knot/link lengths: Fit of the fj states data to tight knot and link 
lengths (ropelengths). Errors are shown for the states, but they are too small to be visible for the 
lengths of ideal knots and links, however we include a 3% error in the knot energies due to the fact 
that we are dealing with physical knots and links. (See discussion in the text.) Non-fitted knots 
and links are not shown. 



We hold the flux $ fixed and vary the field to find the functional form of the energy W 
for a toroidal solenoid. The general form of the energy is 



W = - I B 2 pdzdpda. (1) 



D 



The dzdp integration runs over the cross section of the tube D. Since B is independent of 
a, the a integration gives 



D 



W — 7T J B'pdzdp. (2) 

The flux through D is 

$ = / B dz dp. (3) 



D 



Now we want to vary W with respect to B while holding $ fixed. This is equivalent to 
considering 

5(W - A$) = 0, (4) 

where A is a Lagrange multiplier. For unit vector n normal to the cross section, the variation 
of B gives 

/ {2nBp - An) • SB dz dp = 0, (5) 
Jd 
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FIG. 2. Cross section of a toroidal flux tube of minor radius R\ and major radius R 2 . Both polar 
and cylindrical coordinates are shown. 



which vanishes for arbitrary SB only if 



B{p) 



An 
2wp' 



We find A from the requirement $ = const, which gives 

7T$ 2 



B -- 
W 



with 



dz dp 
d P 



(6) 

(7) 
(8) 

(9) 



To calculate the integral I over the cross section of a torus of major radius R2 and minor 
radius R\, it is convenient to introduce polar coordinates (r, 9) with the origin at the center 
of disk D, plus a toroidal angle a; see Fig. [2j The result of integration over D is 



2tt 



Ri 



which leads to 



W{R 2 



K 2 — ti x 



$2 



i?2 — \j R\ ~ R\ 



(10) 



(11) 



The analogous result for the cylinder of length 2nR 2 is 



R{ 



2 ' 



(12) 
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and so the ratio 



W(R 2 ) {Rx/R 



2 



\2 



W (R 2 ) 2[1 - y/1 ~ (Rl/R2) 2 } 



(13) 



the graph of which is plotted in Fig. [3j will define our new energy. 

Formally, for an embedded tube K of fixed radius R\ and parametric centerline curve 
7(5) with curvature k(s), we define the flux tube energy e(K) by the integral 



<K) = ^IL + J sjl- RW(s) ds J , (14) 

where L = j^ds is the length of the center line. (In numerical studies of tight knots and 
links, it is observed that the integral in (14) is typically ~ which translates into an 



~7% correction of the energy from the ropelength value.) It would be an interesting project 
to numerically minimize the flux tube energy for various knot and link types. However, 
this is likely to be a somewhat challenging project: minimizing functionals of curvature (a 
second derivative of position) constrained by tube contact (a function of position) is quite 
difficult. Still, there has been recent progress in the numerical modeling of elastic rods with 
self-contact |30j which lead to us hope that these computations may be tractable in the near 
future. 

In the meantime, we have chosen to minimize the original Eq(K) energy numerically using 
ridgerunner j27j, and then compute the s(K) energy for the £0 (-^-minimizing configura- 
tions on the grounds that the difference between e(i^)-minimizing and £0 (-^-minimizing 
configurations are likely to be small. In Figure [4] we have histogrammed the shortest 72 
knots and links after curvature corrections have been applied. In Figure [5] we have his- 
togrammed the currently available complete set of 945 curvature corrected knots and links. 

We should remark that in some cases the curvature is discontinuous, but the field is not. 
For example, the inner loop of the chain of three elements (the 2\#2\) is shaped like a 
race track — two straight sections connected by two half circles. In a solenoid of this form 
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FIG. 4. Histogram of the magnitudes of the curvature corrections to the first 72 knot and link 
lengths. 
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FIG. 5. Histogram of the magnitudes of the curvature corrections to all 945 currently tabulated 
knot and link lengths. 



the curvature is discontinuous at the junctions, but if we move along the field lines we find 
that the fields are already changing before they reach the junctions since the windings are 
different in the regions beyond the junction and the field is affected by the fringe fields in 
that region. Hence the fields can be continuous through the junction. 

As an example of a case with both length and curvature correction that can be calculated 
exactly [19J, consider the chain of three unknots 2f ^=2f which has length Qtt + 2 w 20.8496 
and has i?2 = 2Ri in curved regions and R2 — > 00 in straight sections. We find an overall 
corrected value 



4(2 - v/3) 



6tt + 2 « 0.933013(6vr + 2) « 19.4529 



(15) 



Note that all exactly calculable link lengths in our tables are unique, but degeneracies 
can occur at longer lengths. The first such examples are the links corresponding to the Eq 
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and Dq Dynkin diagrams. Both have the length eo(E$) = £q{D§) = 12n + 7 (which falls 
slightly beyond the largest lengths included in our tables) and both have the same curvature 
corrected energy e(E e ) = e(D 6 ) = t^ztm + ^ ■ Other corrections should lift this degeneracy. 



IV. RESULTS 

In our model, the chromoelectric fields [31] are confined to the knotted and linked 
tubes, each carrying one quantum of conserved flux [32j [33] . We consider a stationary 
Lagrangian density 

C=^tiF 0t F 0l -V, (16) 

where, similar to the MIT bag model [34j, we included the possibility of a constant energy 
density V. To account for conservation of the flux we add to £ the term 

trA[$ £ /(vra 2 )-n l F 0l ], (17) 

where n l is the normal vector to a section of the tube of radius a and A is a Lagrange 
multiplier. Varying the full Lagrangian with respect to A^, we find 

D°(F 0l - Xrii) = 0, (18) 
D\F 0t - Xm) = 0, (19) 

which have the constant field solution 

F 0i = (^/vra 2 )^. (20) 

With this solution, the energy is positive and, to first approximation, proportional to 
the length of the tube I and thus the minimum of the energy is achieved by shortening I 
(i.e., tightening the knot), subject to the curvature correction discussed above and other 
corrections discussed below. 

We proceed to identify knotted and linked QCD flux tubes, i.e., curvature corrected 
physical flux tubes, with glueballs and/or predicted glueballs, where we include all / states. 
The lightest candidate is the /o(500), which we identify with the shortest curvature corrected 
knot or link, i.e., the Hopf link 2 2 ; the /o(980) is identified with the next shortest knot or 
link, in this case the trefoil knot 3i, etc. By the fourth knot/link, the ordering begins to be 
reshuffled due to the curvature corrections, see Tables [T] and |III[ 

All knot and link lengths have been calculated for states corresponding to energies well 
beyond 2 GeV. Above ~ 2 GeV the number of knots and links grows rapidly, and so the 
corresponding hadronic states should become dense relative to their typical width. Hence 
we will confine our investigations to knot lengths corresponding to all known / states below 
~ 2 GeV. 

Our detailed results are collected in Table [I] through |IV[ where we list the masses and 
error bars for the / states (other properties can be found in [2]) and our identifications of 
these states with knots and links together with the corresponding knot and link lengths (see 
Tables |l| and [XT]) , curvature corrected lengths (see Tables III, IV and VII) and fitted energies. 



We will give two interpretations of the data. The first possibility is with the /o(1370) 
identified with the 5i knot which results in a prediction of a new state near 1190 MeV 
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FIG. 6. Fit using curvature-corrected lengths: High one parameter /o(1370) fit of the fj states 
data to the curvature-corrected knot and link data. Errors are shown for the states, and the 3% 
estimated from the text is shown for the knots and links. Non- fitted knots and links are not shown 
in this figure. 
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FIG. 7. Predictions using curvature-corrected lengths: The circles are locations of knots and 
links that do not have corresponding fj states in Fig. [6j Hence these are the locations of states 
predicted by the high /o(1370) one parameter fit of the model. 
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FIG. 8. Curvature-corrected lengths: The combined set of measured and fitted states (dots with 
error bars) and predicted states (circles) for the high /o(1370) one parameter fit. 
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FIG. 9. Low /o(1370) one parameter fit with curvature-corrected lengths: This is our best fit of 
the fj states data to the knot and link data. Errors are shown for the states, as is the 3% error 
estimate of knot and link lengths. Non-fitted knots and links are not shown. 
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FIG. 10. Low /o(1370) one parameter fit for curvature-corrected lengths: This figure shows the 
knots and links that do not have corresponding fj states in Fig. [9j Hence these are the locations 
of states predicted by the model. 
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FIG. 11. Low /o(1370) one parameter fit for curvature-corrected lengths: TThis figure shows the 
combined set of measured and fitted states (dots with error bars) and predicted states (circles). 
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identified with the <A\ link. The other possibility, which gives our best fit, is to identify the 
/(1370) itself with the A 2 link to give a one-to-one matching of the first 12 fj states with 
the first 12 knots and links. We will discuss the details of these options below. 

For comparison purposes we have displayed results for knot energies proportional to 
lengths (Tables [T] and [n]) and also for the curvature corrected knot energies, (see Tables 
III, IV and VII). Let us begin with the uncorrected length case. In Figure [I] we display 



a one-parameter least-squares fit to the experimental data (below 1945 MeV) for the mass 
spectrum of / states identified with knots and links. The fit is 

E{K) = A tube e (K), A tubc = 57.14 ± 0.53 MeV, (21) 



where Eq(K) is the dimensionless length of the knot or link K as defined above. The fit (21 ) 
shows fair agreement with our model. One measure of the quality of the fit is given by the 
adjusted R squared, R 2 = .998. Since we have more knots than / states, the relevance of R 2 
must be carefully interpreted and not taken at face value, especially beyond ~ 1700 MeV 
where the lack of particle data per knot becomes pronounced. 

Compared to our 2003 results [TT] . the /o(500) now falls well below the line of our fit and 
we pay a penalty in x 2 . The result is also partially responsible for the change in the value 
of our fit parameter A tu be, compared to what we find using the 2010 values for the /o(600). 
It will be interesting to see if the new /o(500) PDG numbers are stable. Since they are 
extracted from partial wave analysis and multiply subtracted dispersion relations, there are 
questions about comparing with the mass values for other / states found in invariant mass 
plots. There is also the important issue of mixing with four quark and other resonances. 
However, we are not in a position to comment further on these matters with confidence 
except for the general qualitative remark that in a flux tube model we could expect mixing 
with excited mesons if there is resonant behavior of tube breaking (forming a qq pair) and 
rejoining (by annihilating the pair). Resonant breaking and rejoining a two positions, or in 
two link tubes mixes with four quark states, etc. 

Next we consider the curvature-corrected case where the one parameter fit is shown in 
Figure [6j The \ 2 is substantially improved over the uncorrected length case and drops to 
X 2 = 33 to a large extent due to the two most restricting states, the /(1270) and the /(1285), 
have been brought in line with the rest of the fj data by the curvature corrections. Hence 
we continue to consider only the curvature-corrected length case for the remainder of this 
section. 

Figure [7] shows the locations of knots and links with no corresponding / state, hence it 
gives the locations of new states predicted by the model. Our first new state is at 1190 
MeV, corresponding to the <A\ link. It is interesting to consider the PDG entry for the state 
/2(1270). Of the 36 quoted experimental observations, all but one is within three a of the 
PDG average mass of 1275 MeV. The 5.5 a outlier at (1220 ± 10) MeV is from the process 
pp — > ppn + n^ in the experiment of Breakstone et al. [35], and we suggest its identification 
with our predicted state at ~ 1190 MeV. 

We further predict twelve states around (1710 ±20) MeV. To justify such a proposal, one 
only need to look at the PDG entry for the state /o(1710) to see a considerable amount of 
tension in the data with a large number of incompatible mass measurements in this region. 
We interpret this as an indication of multiple J ++ states near 1700 MeV that need to be 
resolved, just as is suggested by our model; see Tables jl] and |III| Similar reasoning applies 
to the states in the vicinity of /2(1810), /o(1910), and ^(1950) where there is also tension 
in the data. A global statistical analysis of /-state data to establish a statistical significance 
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of such an interpretation should be carried out. Figure [8] gives both the combined fit and 
predicted masses and is displayed for convenience. Better HEP data will provide further 
tests of the model and improve the high mass identification. 

We are now ready to discuss the second fit possibility. Since the error on the / (1370) 
mass is rather large (±150 MeV) it can be identified with several different knots and links. 
However, the identifications of the other /-states, except for the /o(500), are much more 
constrained due to the small errors on their masses and this in turn restricts the identification 
of the /o(1370) to two allowed choices with reasonable x 2 s. We call the case discussed above 
the high-fit, where we identify the /o(1370) with the 5i knot, and the low-fit where we 
identify the / (1370) with the 4 2 link. The high-fit is the one that predicts a new state at 
~ 1190 MeV, while the low-fit gives a one-to-one match between the first 12 /-states and 
the first 12 knots and links. The low-fit gives a somewhat better x 2 , but the high-fit gives an 
acceptable \ 2 an d could be required if for instance the 1220 MeV state of [55] is confirmed 
making it necessary to free up a low energy knot/link to identify with it. Since we have 
already presented the high-fit above, we now proceed to discuss the low-fit. 

The low-fit leaves no gaps in the spectrum until we get near 1700 MeV. The x 2 is improved 
and the fitted value of A tube is similar to the high-fit case. We find 

E(K) = A tuhc e(K), A tube = 65.16 ± 0.61 MeV. (22) 



The low-fit (22) is our best overall fit to the data. We have tabulated the fitted energies for 
the first 72 knots and links, identifying all /-states of energy less than 2 GeV, hence giving 
predictions of many new states at 1690 MeV and above. 

Now let us discuss the statistical significance of our results. First we have done a number 
of tests to determine if our hypothesis that glueballs are knotted flux tubes is likely to be 
correct. Let us begin with the well-known Kolmogorov-Smirnov test, which tests distribution 
functions F(y) and F(y) by means of the quantity 

axtp\F(y)-F(y)\, 
y 

which in this case are the distribution functions of the glueball masses and calculated knot 
and link lengths. The resulting p-value for the Kolmogorov-Smirnov test is 

Pks = 0.95, (23) 

where we recall that p is bounded < p < 1 and p < 0.01 implies poor correlation, 
0.01 < p < 0.05 implies moderate correlation and 0.1 < p implies strong correlation. Hence 
the Kolmogorov-Smirnov test implies our model is in excellent agreement with the data. 

The Kolmogorov-Smirnov test is a measure of goodness-of-fit. We summarize this and 
a number of other goodness-of-fit tests as well as several variance tests in Table |Vj All 
show excellent agreement between model and data. Here and below we give p- values for the 
high-fit case. The low-fit values are similar. 

Another approach is to calculate the x 2 f° r the data set, subject to the corrections of 
energy per unit length differences and deformations from the ideal knot case. We have 
argued that the minimum energy and minimum length of knots do not necessarily coincide. 
While we expect the average energy per unit length of a knot does not strongly depend on 
knot type, there is still a small correction due to this effect. In addition, the tubes can be 
constricted due to being wrapped by another section of the tube or distorted by wrapping 
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FIG. 12. Cross section of a toroidal flux tube of minor radius R± and major radius R 2 constricting 
a cylindrical flux tube of radius R(z). 



tightly around another section of tube (as a rope wrapped tightly around a post). We can 
approximate such corrections and will consider a typical example below. 

We have calculated the change in energy of a pair of linked flux tubes to provide an 
example of corrections we should expect due to constriction, which in turn gives a contri- 
bution to the expected error in using physical tubes instead of ideal mathematical tubes to 
model glueballs. The example we consider is a cylindrical tube along the z-axis encircled by 
a toroidal tube lying in the xy-pla.ne at z — 0. We assume both tubes have circular cross 
section with that of the torus staying fixed, but that of the cylinder being of reduced radius 
in the region of constriction, 



R(z) = R(0) +Ri- (Rf - z 2 ) 



2\l/2 



where 



Zrn ^ Z ^ Zin 



[R\ - R 



1/2. 



(24) 



(25) 



see Fig [12] We further assume both tubes carry the same amount of flux and so set their 
undistorted radii equal, R\ = R(z m ). The torus would tighten until R 2 = R±, except that it 
begins to encounter the cylindrical tube at R2 = 2R\. The torus energy 



AW 1 = W{R 2 ) - W(2R 



(26) 



where the function W is defined in ([3]), falls as it tightens and the energy in the cylindrical 
tube 



AW 2 



$1 



7T 



Zm 
K 2 



dz 



(27) 



, R?{z) 

grows as it is constricted. Stability is reached by minimizing AW\ + AW 2 with respect to 



R(0); see Fig. 12 



We have estimated the constriction for the Hopf link of two magnetic flux tubes and find 
a ~ 30% correction over the region of constriction, which translates into an overall ~ 5% 
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correction to the link energy, since about 15% of the Hopf link is constricted. We expect 
distortion effects to be similar. However, since all knots and links have similar corrections 
that modify their total energies in the same direction, we expect the spread in variation to 
be smaller than the correction itself. Hence a 5 ~ 5% error on physical knot energies versus 
tight knot energies is not unreasonable. An actual QCD flux tube could be more rigid due 
to confinement effects, so we can justify reducing the error bars to something smaller. Even 
assuming QCD flux tubes are substantially more rigid than magnetic flux tubes, we still 
find an acceptable value for x 2 . For example, letting 5qcd = 3% we find the model is in 
reasonable agreement with the data as seen in the figures for the fits. 

In terms of the bag model [3lj, the interiors of flux tubes of tight knots correspond to the 
interiors of bags. The flux in the tube is supported by current sheets on the bag boundary 
(surface of the tube). Knot complexity can be reduced (or increased) by unknotting (knot- 
ting) operations [361 El] ■ In terms of flux tubes, these moves are equivalent to reconnection 
events [38]. Hence, a metastable glueball may decay via reconnection. Once all topologi- 
cal charge is lost, metastability is lost, and the decay proceeds to completion. Two other 
glueball decay processes are: flux tube or string breaking JSHHH] (this favors large decay 
widths for configurations with long flux tube components) and quantum fluctuations that 
unlink flux tubes (this would tend to broaden states with short flux tube components). Since 
the publication of some minor quantitative progress has been made in understanding 
knot flux tube decay, but these results are still insufficient to go beyond the qualitative 
observations made in [TT] . 

Let us make one final comment about the model. A more conservative approach is to 
assume only the ++ states, i.e., the f states, correspond to knotted/linked flux tubes. 



In that case there are only five states to fit and the result is displayed in Figure 13 for 
comparison, where the identified states are 

[/o(500),/o(980),/o(1370),/ (1500),/o(1710)] <— > [2 2 , 3 X , 5 X , 5 2 , 6 2 ]. (28) 

The slope is essentially unchanged, the R 2 = .998 value is roughly the same. The masses of 
the predicted states can be easily gotten by rescaling the knot lengths in the table with the 
new slope parameter A tubo = (65.50 ± 1.81) MeV. Note that the / (1370) has been moved 
in the ordering to improve the fit. More data is needed to distinguish between the fit of all 
fj states and the restricted fo fit. 

The x 2 = 15.62 for the fit is not particularly good but it only takes replacing the new 
/o(500) values with the old / (500) mass and error bars to get \ 2 — 0.56. So if the new 
/o(500) mass is due to mixing and one could extract the unmixed value for the pure gluonic 
state it is possible that the fit would improve again. 

If a sufficient number of fj states are found, so that they outnumber the total number 
of knots and links, then this would be evidence to support the restricted fit (since all short 
knots and links are presumed to be known). However, this is not the case at present. 



V. DISCUSSIONS AND CONCLUSIONS 

We have considered hadronic collisions that produce some number of baryons and mesons 
plus a gluonic state in the form of a closed QCD flux tube (or a set of tubes). From an 
initial state, the fields in the flux tubes quickly relax to an equilibrium configuration, which 
is topologically equivalent to the initial state. (We assume topological quantum numbers 
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FIG. 13. The /o states data is fitted to the curvature corrected knot and link data. Errors are 
shown for the states and estimated to be 3% for the knot/link energy. Non-fitted knots and links 
are not shown. 
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FIG. 14. The fo states data is fitted to the curvature corrected knot and link data and predicted 
states (open disks). 
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are conserved during this rapid process.) The tube radius is set by the confinement scale, 
so to lowest order the energy of the final state depends only on the topology of the initial 
state and equals the length l K of the tube times the average energy per unit length, or 
the dimensionless knot or link length Eq(K) times the energy scale parameter A tu b e - While 
related to Aq CD by constants of order unity, A tube can be more accurately determined (see 
above) and hence could be a useful dimensionful parameter in studying other properties 
of QCD, such as scattering and hadronization processes. The relaxation proceeds through 
minimization of the field energy. This process occurs via shrinking the tube length and the 
process halts to form a tight knot or link. Flux conservation and energy minimization also 
force the fields to be homogeneous across the tube cross sections for straight tube sections 
and the fields fall like 1/p for curved sections as shown above. We have estimated corrections 
to the simple energy-length proportionality and have used them to correct and place error 
bars on the physical knot lengths. 

Details of knot excitations would be interesting to investigate, as would other quantum 
corrections, but at present we do not have a reliable way to estimate these effects, nor do 
we have a good way to calculate glueball decay rates. However, we do expect high mass 
glueball production to be suppressed because more complicated non-trivial topological field 
configurations are statistically disfavored and we also expect higher mass glueballs to be 
relatively less stable. 

On the lattice the glueball is associated with a plaquette operator that, when operating 
on the vacuum, creates a closed loop of chromoelectric flux [12]. This loop is a path on 
a square lattice and the glueball mass should be proportional to the length of the path. 
Using our assumption that topologically trivial paths are too unstable to allow measurable 
masses and assuming that we are studying single flux tubes on the lattice (no links), the 
first stable closed loop on the lattice will be the trefoil. The shortest length for a trefoil on a 
square lattice is 24 lattice spacings (24 tube diameters in dimensionless units), so, without 
smoothing, the lattice should predict the lightest glueball mass to be a factor of 24/e(3i) 
larger than the value we quote in the table, i.e., 

24 

£iattice(3i) ~ ^tt£(3i) ~ 1450 MeV. (29) 
e{6i) 

While this naive result is only a rough approximation, it could be refined and does indicate 
that lattice calculations of glueball masses can come out on the high side. It is also an 
explanation of why our results differ from lattice predictions. 

In addition to not fitting naturally into the quark model [43J, glueballs have some other 
characteristic signatures, including enhanced production via gluon rich channels in the cen- 
tral rapidity region, branching fractions incompatible with qq decay, very weak coupling to 
77, and OZI suppression. All the /-states have some or all of these properties. For instance, 
none have substantial branching fractions to 77. However, mixing with qq isoscalar states 
can obscure some of these properties. All these observations are in qualitative agreement 
with the model presented here. 

Our high-fit model predicts one new state at 1190 MeV, twelve states concentrated near 
1700 MeV and a tower of higher mass states with the next dense concentration starting near 
1900 MeV. The low-fit model makes similar prediction except that there is no new state 
near 1200 MeV. We have argued that there is sufficient tension in the experimental data in 
these regions to allow the identification of many more states with knots and links. A careful 
statistical analysis of the data of all /-regions to resolve hidden states is needed. Recall we 
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FIG. 15. Uncorrected lengths: Two parameter fit of the fj states data to the knot and link 
data. Errors are shown for the states, but they are too small to be visible for the knots and links. 
However, a 3% knot length error is included for reasons discussed in the text. Non-fitted knots 
and links are not shown. 



are assuming that J is intrinsic angular momentum and not rotational angular momentum 
as we assumed in [TT] . 

As a variant example of the models we have been considering and as a comparison, we 
consider a two-parameter fit to the fj data where the origin is not fixed at zero glueball 
mass and zero tube length. Fitting the non-curvature-corrected length data to the glueball 



data gives a shallower slope and an intercept at positive glueball mass; see Fig. [15} While 
the x 2 is somewhat better than for the non-curvature-corrected one-parameter fit, we do 
not have an interpretation of the non-zero intercept other than a zero length tubular bag 
constant which seems rather unphysical. This problem does not arise for either the high- or 



low-fit two-parameter models with curvature-corrected lengths. As seen in Figs. [16] and [T7 
the intercepts in both these models are consistent with zero. This result and in addition to 
the improved \ 2 when we use the curvature-corrected knot lengths in the one-parameter fits 
gives us confidence that either our one-parameter (i.e., the slope A tu & e ) high or one-parameter 
low curvature-corrected fits are the sufficient and best choices for a robust model. 

We note that some of the states can be degenerate at lowest order due to the symmetries of 
some knots and links. The trefoil comes in two versions, left- and right-handed. Some prime 
knots have aZ 2 xZ 2 symmetry leading to a potential four-fold degeneracy. The cases of links 
and composites are even more complicated, with a variety of higher degeneracies possible 
[H]. How these degeneracies can be lifted at higher order is a topic for future research, but 
let us make a few comments here. Since many knots and links come in more than one form, 
we can consider their mixings. For example, the trefoil can be left-handed 3^ or right-handed 
3f . Since the trefoil can decay by changing a single crossing, changing two crossings can take 
3^ to 3f . We have identified the 3i with / (980) whose width is ~ 50 MeV. Assuming the 



20 




s(K) 



5 10 15 20 25 30 35 

FIG. 16. Two parameter fit with curvature-corrected length. Two parameter low /o(1370) fit of 
the fj states data to the knot and link data. Errors are shown for the / states, and a 3% error is 
estimated for knot and link lengths. Non-fitted knots and links are not shown. 




e(K) 



5 10 15 20 25 30 35 

FIG. 17. Two parameter fit with curvature-corrected length. Two parameter high /o(1370) fit of 
the fj states data to the knot and link data. Errors are shown for the / states and estimated for 
the knots and links. Non-fitted knots and links are not shown. 
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decay is via reconnection or quantum tunneling, then we expect the 3^ and 3f versions to 
mix at the 10~ 3 level by changing two crossings. This corresponds to a ~1 MeV level splitting 
and is undetected experimentally. As the quantum state = A= (3\ — 3f) has self-linking 

number SL(^ 3l ) = 1 when properly normalized, while if)' 3l = ^ (3^ + 3f) has SL(?/> 31 ) = 0, 
we would expect the tpL to be the more stable linear combination, and the proper state to 
be identified with the /o(980). Similar comments apply to all our identifications. Note that 
the Hopf link comes in only one form, so it corresponds to the single broad state /o(500). 
Mixing could potentially be detected for other broad states where splitting could be a few 
MeV, but only when they are identified with knots or links that come in more than one 
type. 

It would be trivial to extend our predictions to states above 2 GeV. One just takes the 
knot of appropriate lengths from [27] and scales them by the dimensionful parameter A tu be 
from the fit. However, there is insufficient mesonic data above 2 GeV to improve or constrain 
our fit. More experimental data to test the model in this region would be very welcome. 

Knotted and linked magnetic fields configurations have been discussed with respect to 
a number of plasma phenomena including spheromaks |45j, astrophysical [46J and atmo- 
spheric jUJ objects. Similar comments apply to Bose-Einstein condensates |48j . and various 
field theories. The system we consider is intrinsically different since we assume that in QCD 
we have confinement into tubes and then the tubes get knotted and linked, as opposed to 
finding knotted or linked fields that then may or may not get confined. In our case we 
expect energy proportional to tube length (with corrections as discussed in this paper), 
where the later case not involving initial confinement would not necessarily be expected to 
have a simple length-energy relationship. The energies and sizes of these classical solitons 
are sometimes difficult to quantify since they depend on parameters of the plasma, includ- 
ing temperature, pressure, density, ionic content, etc. However, it was argued in [TTJ that 
well-defined topological soliton energies can be identified for vacuum QCD or any vacuum 
quantum-flux-tube system. Hence we emphasize that in all systems which support quantum 
flux tubes, including those occuring in media, from the quark-gluon plasma to superconduc- 
tors, the energy spectrum of knot and link solitons will be universal up to a scaling for fixed 
values of parameters (IHl [50] . 

To conclude, we have given two interpretations of the /-state data with a model of 
knotted chromoelectric flux tubes in QCD. The first possibility is where the /(1370) is 
identified with the 5i knot which results in a prediction of a new state around 1190 MeV 
which is identified with the A\ link. The other possibility — which also gives our best fit — is 
to identify the /(1370) with the A\ link to give a one-to-one matching of all the first twelve 
fj states with the first twelve knots and links. Experiments could help to resolve which of 
these two possibilities is the correct choice. 

Finally we should point out that there is considerable amount of tension in the fj data, 
as indicate by the x 2 s of the individual states quoted by the PDG. We would not expect our 
fits to be better than the fits of the data on which they are based. 
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TABLE I. Comparison of the glueball mass spectrum and fit energies for £q(K) less than ~ 36 
E\(K) and E2{K) are for one and two parameter high /q(1370) fits. 
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24.893 


1422 


1421 


A (1510) 


1518 ±5 


6 


', 

3 




25.181 


1439 


1434 


/a (1525) 


1525 ±5 


6 


2 

1 




27.146 


1551 


1525 


A (1565) 


1562 ± 13 


7 


2 
7 




27.760 


1586 


1554 


A (1640) 


1639 ± 6 


(2 T #2f#22) kc 


28.133^ 


1608 


1571 






2?#4f 


28.311 


1618 


1579 






6 


2 
2 




28.356 


1620 


1582 






6i 




28.364 


1621 


1582 






62 




28.522 


1630 


1589 






3m#3^j 


28.537 


1631 


1590 






2? #4i 


28.742 


1642 


1599 






7 


2 
S 




28.886 


1651 


1606 






6 


3 
1 




28.914 


1652 


1607 






63 




28.929 


1653 


1608 






6^ 




29.006 


1657 


1612 






6 


2 




29.057 


1660 


1614 






(2?#2?#2?) lc 


29.133^ 


1665 


1618 


[/o(1710)]i 


[1720 ± 6] i 


8 


3 
7 




30.297 


1731 


1672 






819 




30.502 


1743 


1681 






7i 




30.715 


1755 


1691 


[A(1710)] 2 


[1720 ±6] 2 


820 




31.557 


1803 


1730 


[A (1810)] i 


[1815 ± 12] i 


2?#5i 


31.908 


1823 


1746 






7 3 




31.975 


1827 


1749 






8? 5 




32.093 


1834 


1755 






72 




32.122 


1836 


1756 






7 


2 
1 




32.129 


1836 


1756 






7 4 




32.146 


1837 


1757 






3i m #4f 


32.189 


1839 


1759 






8 


3 
3 




32.514 


1858 


1774 






7-> 
' 2 




32.520 


1858 


1775 






7-' 

'4 




32.542 


1860 


1776 
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TABLE II. Comparison of the glueball mass spectrum and fit energies for e{K) less than ~ 36. 



State 



Mass 



A 



e(K 



E 1 (K 



E 2 {K 



[/ 2 (1910)]i,[/ 2 (1810)] 2 
[/2(1950)]i 



[1903 ±9]i, [1815 ± 12] 2 
[1944 ± 12] i 



3i#4i 
7 5 

7 2 
'3 

°10 
8 21 

7 7 
7 6 

n 

7? 

2 2 #2 2 #3 1B 
949 

2?#5 2 
9« 



8 i(i 
Q 2 

8? 



946 

Q 2 

Q 2 

a 61 
942 
Q 2 



[/ 2 (1910)] 2 



[1903 ± 9] 2 



9 2 



(2 2 #2 2 #2?#2 2 ) fcc 



9 2 
•'r>4 

10124 



9 44 
9 56 



[/2(1950)]s 



[1944 ± 12] 2 



32.636 
32.641 
32.659 
32.675 
32.735 
32.771 
32.816 
32.857 
32.873 
32.911 
32.959 
33.028 
33.041 
33.134 
33.167 
33.221 
33.359 
33.711 
34.008 
34.214 
34.319 
34.350 
34.689 
34.755 
35.083 
35.276 
35.416 
35.516 
35.526 
35.555 
35.745 
35.844 
35.895 
36.080 



1865 
1865 
1866 
1867 
1871 
1873 
1875 
1878 
1878 
1881 
1883 
1887 
1888 
1893 
1895 
1898 
1906 
1926 
1943 
1955 
1961 
1963 
1982 
1986 
2005 
2016 
2028 
2029 
2030 
2032 
2043 
2048 
2051 
2062 



1780 
1780 
1781 
1782 
1785 
1786 
1788 
1790 
1791 
1793 
1795 
1798 
1799 
1803 
1805 
1807 
1813 
1830 
1844 
1853 
1858 
1859 
1875 
1878 
1893 
1902 
1912 
1913 
1914 
1915 
1924 
1929 
1931 
1940 



Notation n k means a link of I components with n crossings of k type, see e.g 



K#K' stands for 



the connected sum of K and K' , and (K) m is the mirror image of K. 

All values of the ropelength sq(K) are from [27], except for composite links that were calculated 
separately but also with Ridgerunner. The dimensionless length convention agrees with jllj and is a 
factor of 2 smaller than that in [57] . 

E(K) is the value of the fitted mass corresponding to the PDG or predicted mass obtained from e(K). 
Dimensionless curvature corrected knot energy. 

Exact values: sq(2\ ) = 47r (Hopf link), e (2f #2f ) = 6tt + 2 (chain of three links), 



e ((2?#2f #2?) kc ) = 8tt + 3 (key chain link with three keys), £ ((2?#2f #2f ) 



B7r + 4 (linear chain 



with four links), and e ((2f #2f #2f #2f ) kc ) = 10n + 4 (key chain link with four keys). 
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TABLE III. Low /o(1370) curvature corrected fit: Comparison of the glueball mass spectrum and 
fit energies for e(K) less than ~ 32. 



State 



Mass 



e(K' b 



/o(500) 
/o(980) 
/o(1370) 
/ 2 (1270) 
/i(1285) 
/l(1420) 
/ 2 (1430) 
/o(1500) 
/i(1510) 
/a (1525) 
/ 2 (1565) 
/ 2 (1640) 



/o(1710) 



/ 2 (1810) 



475 ± 75 
990 ± 20 
1350 ± 150 
1275.1 ± 1.2 
1282.1 ± 0.6 
1426.4 ± 0.9 
« 1430 
1505 ± 6 
1518 ± 5 
1525 ± 5 
1562 ± 13 
1639 ± 6 



1720 ± 6 



1815 ± 12 



2 2 
3i 

1? 

4i 
2?#2? 

5i 
2?#3i 

5 2 

5? 

6§ 

6? 

7? 

61 

6i 
2?#4f 

3l#3l m 

3i#3i 
62 
6? 

(22#22#2f) kc 
2?#4x 
63 
65 
7i 

6 :i , 

(22# 2 2#2f) lc 
7i 

8? 
Si!) 

7 3 
820 
7? 
7 2 
2?#5i 

7 4 
3lm#4? 

8? 5 
3i#4? 

7 2 
'4 



11.724 
14.943 
18.250 
19.411 
19.556 
21.559 
22.697 
22.779 
22.866 
23.309 
24.854 
25.735 
25.924 
26.025 
26.046 
26.135 
26.151 
26.158 
26.327 
26.449 
26.466 
26.567 
26.590 
26.720 
26.963 
27.449 
28.018 
28.152 
28.458 
29.025 
29.151 
29.231 
29.330 
29.339 
29.385 
29.402 
29.496 
29.536 
29.544 



764 
974 
1189 
1265 
1274 
1405 
1479 
1484 
1490 
1519 
1619 
1677 
1689 
1696 
1697 
1703 
1704 
1704 
1715 
1723 
1724 
1731 
1733 
1741 
1757 
1788 
1826 
1834 
1854 
1891 
1899 
1905 
1911 
1912 
1915 
1916 
1922 
1924 
1925 
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TABLE IV. Comparison of the glueball mass spectrum and fit energies for e{K) less than ~ 32. 



















State 


Mass 


K 


a 


e(K\ 


b 


E X (K\ 


C 



22#2f#3iA 29.682 1934 

3i#4i 29.790 1941 

7 5 29.806 1942 
/ 2 (1950) 1944 ±12 2j#5 2 29.840 1944 

7§ 29.873 1946 

7 6 29.894 1948 
7| 29.895 1948 

22#5? m 29.929 1950 

2f#5f A 29.952 1952 

8| 29.957 1952 

7§ 30.015 1956 

821 30.017 1956 

7 7 30.092 1961 
7? 30.112 1962 
7§ 30.302 1974 
9| 3 30.416 1982 
8f 6 30.525 1989 
8| 30.571 1992 
8f 30.605 1994 

2f#2j#3i S 30.611 1995 

9| 9 30.839 2009 

8^ 30.967 2018 

8f 31.214 2034 

8| 31.473 2051 

9 4 6 31.513 2053 

9§ 31.521 2054 

9;?3 31.776 2070 

9§ x 31.909 2079 

9| 7 31.916 2080 

9 42 31.950 2082 



a Notation n' fe means a link of ^ components with n crossings of fc th type, see e.g., [36]. K#K' stands for 

the connected sum of K and K' , and (K) m is the mirror image of K. 
b All values of the curvature corrected knot energies 4>(K) are modified from [27], except for composite 

links that were calculated separately but also with Ridgerunner. The dimensionless length convention 

agrees with [TT] and is a factor of 2 smaller than that in [27] . 
c Ei{K) is the value of the fitted mass corresponding to the PDG or predicted mass obtained from e(K). 
d Dimensionless curvature corrected knot energy. 

e Exact values: £ (2?) = 4vr (Hopf link), e (2?#2?) = 6tt + 2 (chain of three links), 
e ((2?#2?#2?) kc ) = 8tt + 3 (key chain link with three keys), e ((2?#2?#2?) lc ) = 8tt + 4 (linear chain 
with four links), and e ((2f #2f #2f #2f ) kc ) = 10tt + 4 (key chain link with four keys). 
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TABLE V. Statistical tests of the model. Recalling that p is bounded < p < 1 and p < 0.01 

implies poor correlation, 0.01 < p < 0.05 implies moderate correlation and 0.1 < p implies strong 
correlation, we see that all these tests strongly support the model. 

Goodness-of-fit test p-value Variance test p-value 

Pearson \ 2 0.66 Brown-Forsythc 0.74 

Kolmogorov-Smirnov 0.95 Fisher Ratio 0.69 

Cramcr-von Mises 0.96 Levene 0.74 

Anderson-Darling 0.97 Siegel-Tukey 0.82 
Kuiper 0.99 
Watson U Square 0.90 



TABLE VI. Fit parameters for the model at 95% CL. We collect the results for the various choices 
of one and to parameter fits so raw length fits can be compared with curvature corrected fits. 

Fitting Parameter lp-lcngth 2p-length lp-cc-high 2p-cc-high lp-cc-low 2p-cc-low 

X 2 83.9 45.4 33.3 33.1 28.3 28.4 

Atube 57.14 ± 0.53 46.36 ± 2.64 65.06 ± 0.61 63.83 ± 3.57 65.16 ± 0.61 64.34 ± 3.59 

Intercept 267.11 ± 69.55 28.74 ± 82.18 19.16 ± 82.51 

R 2 0.998 .967 .999 .977 .999 .980 



TABLE VII. High /o(1370) curvature corrected fit: Comparison of the glueball mass spectrum 
and fit energies for e(K) less than ~ 32. Except for the /(4^), this table contains only the fitted 
PDG states. Predictions for other states can be gotten by multiplying the knot energy by the 
appropriate value of A tu b e - 



State 


Mass 


K 


e(K) 


Ei(K) 


/o(500) 


475 ± 75 


2? 


11.724 


763 


/o(980) 


990 ± 20 


3i 


14.943 


972 


/(4?) 




4? 


18.250 


1187 


/ 2 (1270) 


1275.1 ± 1.2 


4i 


19.411 


1263 


/i(1285) 


1282.1 ± 0.6 


2?#2 2 


19.556 


1272 


/i(1420) 


1426.4 ± 0.9 


5i 


21.559 


1403 


/o(1500) 


1505 ± 6 


2 2 #3i 


22.697 


1477 


/ 2 (1430) 


sj 1430 


5 2 


22.779 


1482 


/i(1510) 


1518 ±5 


5? 


22.866 


1488 


#(1525) 


1525 ± 5 


6§ 


23.309 


1516 


/ 2 (1565) 


1562 ± 13 


6? 


24.854 


1617 


/ 2 (1640) 


1639 ± 6 


7? 


25.735 


1674 


/o(1370) 


1350 ± 150 


61 


25.924 


1687 


/o(1710) 


1720 ± 6 


(2 2 #2 2 #2 2 ) kc 


26.449 


1721 


/ 2 (1810) 


1815 ± 12 


7i 


28.018 


1823 


/ 2 (1910) 


1903 ± 9 


7 2 


29.330 


1908 


/ 2 (1950) 


1944 ± 12 


2 2 #5 2 


29.840 


1941 
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